Phase Dynamics of Entangled Qubits 
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We make a geometric study of the phases acquired by a general pure bipartite two level system 
after a cyclic unitary evolution. The geometric representation of the two particle Hilbert space 
makes use of Hopf fibrations. It allows for a simple description of the dynamics of the entangled 
state's phase during the whole evolution. The global phase after a cyclic evolution is always an entire 
multiple of ir for all bipartite states, a result that does not depend on the degree of entanglement. 
There are three different types of phases combining themselves so as to result in the nn global phase. 
They can be identified as dynamical, geometrical and topological. Each one of them can be easily 
identified using the presented geometric description. The interplay between them depends on the 
initial state and on its trajectory and the results obtained are shown to be in connection to those 
on mixed states phases. 
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I. INTRODUCTION 

Quantum phases are a peculiarity of quantum physics 
giving it some of its intriguing and unusual effects, such 
as interference. While relative phases are extremely im- 
portant for any measurable property of a given quantum 
system, global phases are, in general, irrelevant and play 
no major part. Nevertheless, they can be measured using 
conditional dynamics. Supose one has an initial state de- 
scribed by 1 /a/2 I \&) x (|0) 2 + 1 1)2), where the subscripts 
refer to two different subspaces. Suppose now that, if 
particle 2 is in state |1) 2 , \$>) 1 gains a phase of <f>. By 
rotating the second particle so as |1) — > -4=(|1) + |0)) 



and |0) 



f (|0) — |1)), we have that our state becomes 



1 

!)(!*)! |0) 2 + (e* - !)(!*>! |1) 2 ). Detecting the 



probability of finding the second particle in state |1) 2 
gives us P2 oc 1 — cos <p allowing for the determination of 
<f>. 

Based on this discussion, we can now define the total 
phase gained by a system after its evolution. The defini- 
tion used here is the same as proposed by Pancharatnam 
PJ. In his seminal work, he defined the phase of a quan- 
tum state |<&) relative to another state \^f) as 



arg(*|$) 



(1) 



This definition is called from now on the total (global) 
phase. The total global phase gain may originate from a 
combination of phase effects of different types. One possi- 
ble origin of the global phase gain is of dynamical nature. 
This phase is given by the eigenvalues of the Hamiltonian, 
determining the time evolution of its eigenstates. How- 
ever, dynamical phases are not the only type of phase 
appearing during a state's evolution. In order to be more 
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precise, let's restrict ourselves to a two level system, or 
a qubit. The Hilbert space of such a system has a cer- 
tain geometry that can be described by the Bloch sphere. 
This geometry accounts for geometric phase gains. Such 
phases were first introduced in quantum mechanical sys- 
tems by M. Berry 00 and can lead to some amazing and 
counter intuitive effects associated to particles and their 
trajectories. One of the most famous of such effects is the 
Aharanov-Bohm effect, in which a measurable meaning 
is given to the vector potential. The geometric phase de- 
pends on the trajectory realized by the qubit under the 
action of an evolution operator. For cyclic evolutions, it 
is given by half the solid angle of the closed path traced 
in the Bloch sphere, i.e., the area of the surface enclosed 
by the trajectory. 

The development of quantum mechanics motivated 
new questions and problems concerning geometric 
phases. These problems involve the study of geometric 
phases in many particle entangled systems such as 
twin photons @ and Bose-Einstein condensates [(|. Lo- 
cal and non local aspects of geometric phase have also 
been studied 0, as well as its behavior in dissipativc 
systems More recently, it was found that geometric 
phases can also be useful in quantum information theory. 
It was shown that conditional quantum logics, a funda- 
mental part to realize universal quantum computation, 
can be done using geometric phases a result already 
demonstrated experimentally • The interest of these 
so called topological gates is that they are more robust 
and naturally fault tolerant [111 ]. 

These recent results together with the fundamental in- 
terest of the subject, are a motivation for the present 
work: we study here the geometric representation of the 
phase dynamics of pure two-level bipartite systems with 
an arbitrary degree of entanglement. For this kind of 
system, a third possible type of global phase gain can 
be identified: a topological phase, which is a consequence 
of the geometry of the entangled two-level system. This 
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phase has been studied in |12( for maximally entangled 
states (MES) and it is at the origin of singularities ap- 
pearing in the phase of MES during a cyclic evolution y[ ■ 
The topological origin of this phase can be made clear by 
applying the geometric representation of pure entangled 
qubits presented in [Tflj . As shown in the present paper, 
the use of the same representation allows for a clear geo- 
metric picture showing the origin and interplay between 
the three mentioned different types of phases during the 
state's evolution. It is shown that, after a cyclic evo- 
lution, the combination of the phases presented above 
always lead to a global phase gain of an entire multiple 
of 7T. This result, already known and verified experimen- 
tally 01 for a single qubit is recovered here for entangled 
qubits with an arbitrary degree of entanglement. 

The paper is organized as follows: we start by studying 
the case of one qubit and analogous systems, investigat- 
ing their geometric representation and phases acquired 
during a cyclic evolution. We do so to set notations and 
notions for the geometric representation and phase dy- 
namics of two qubit systems, presented in Section III. 
Finally, Section IV is devoted to the study of some ex- 
amples of evolutions leading to different types of phase 
gains for entangled systems. We finish the paper with 
some concluding remarks and a brief discussion on how 
the presented ideas could be tested experimentally. 



II. ONE QUBIT: GEOMETRY AND PHASES 
A. Pure states 




1°) 

FIG. 1: Bloch sphere representation of an one qubit state. 

A general pure one qubit system can be represented 
by a S 3 sphere in i? 4 . S 3 Hopf fibrations define a map 

S 3 — ► S 2 , "decomposing" the S 3 sphere into a S* 2 base 
and a fiber. It has been shown [l3L ITEj that this way of 
representing the space of one qubit is analogous to the 
Bloch sphere: the base can be identified to the Bloch 



sphere while the fiber represents the global phase inde- 
termination of such a representation. Taking the general 
pure state |\&(0)) = a |0) + /3 11), where a and (3 are com- 
plex coefficients, its representation in the Bloch sphere 
is shown in Fig. QJ. The state's coordinates can be 
obtained as follows: 

X = (a x )y = 2Re(a*P) 

Y = (<t v )« = 2Im(a*/3) (2) 

Z=(* t )v = M 2 - |/3| 2 . 

(3) 

As a consequence of the normalization condition, the 
radius of the Bloch sphere is equal to one. Rotations 
around axes in an arbitrary direction displace the state 
vector on the Bloch sphere. We can study, for example, 
the effect of a rotation around the z axis, produced by 
the Hamiltonian H = hwa z /2. Physically, this may cor- 
respond to a magnetic field applied along the z axis in 
the case where our qubits are spin 1/2 particles. At time 
t, the state reads \^(t)) — e^ 1 a\0) + /3 1 1} , where 
u> is the frequency of the applied field. In the Bloch 
sphere, this means that the state actually makes a pre- 
cession around the z axis, reaching again the initial state 
for t = 2n/(J. This evolution leads to a phase gain. To 
compute the global phase gain of the state at any time t 
and any point of the trajectory, one can use Pancharat- 
nam's definition (JTJ. We see that after the 2ir rotation 
the global phase gain is equal to 7r, irrespectively of the 
choice of the initial state. Nevertheless, the nature of this 
phase depends on the initial state and on the Hamilto- 
nian acting on it. Taking, for example, \a\ = cos (8/2) 
and \P\ = sin (0/2), one can identify, in this simple one 
qubit dynamical phase 4>d, defined as 

4>d= I (U f U)dt, (4) 
Jo 

where U is the evolution operator associated to the 
Hamiltonian H . For the particular example studied, we 
have <j>d = — 7rcos6>. As can be seen, there is a depen- 
dence on the angle 9 that characterizes our initial state. 
We can also see that the dynamical phase does not ac- 
count for the total phase gain. In order to get a complete 
description of the total phase gain, one should also con- 
sider the geometric phase <fi g , which is a property of the 
space of states. In the case of qubits describing trajecto- 
ries in the Bloch sphere, the cyclic geometric phase cor- 
responds to the area enclosed by the trajectory realized 
by the state vector on this same sphere. It assumes the 
value of 4> g = — 7r(l — cos#) in our simple example. As 
can be easily checked, we have <f>t = (j)g + 4>d = ~ ?r after 
a cyclic evolution, and this result is state independent. 
The initial state determines only the proportion of the 
geometrical and dynamical phases. 
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B. Mixed states 

One qubit mixed states are described by a density ma- 
trix of the form 

p = a Q0 |0) (0| + ooi |0) (1| + a* w |1) (0| + a n |1) (1| (5) 

They can be geometrically represented in the Bloch ball, 
a generalization of the Bloch sphere. The Bloch ball has 
a radius smaller than one. It is, in fact, a function of 
the purity P = 1 — Trp 2 of the system, as follows: r = 
|1 — 2P] 1 / 2 . P assumes the value of 1/2 for completely 
mixed states and 1 for pure states. The coordinates of 
(0) on the Bloch ball are given by 

X = {a x ) 9 = 2Re(a 01 ) 

Y = (o- y )» = 2Im(a i) (6) 

Z = (cr-)* = °oo - ill- 
It is important to notice that the representation of a 
mixed state in the Bloch ball is different from the pure 
case. A point in the Bloch ball determines the density 
matrix with no global phase ambiguity. Also, the phases 
involved in a cyclic evolution of a mixed state deserves 
more attention. The global phase, analogous to the Pan- 
charatnam's phase, is given by (ftf 1 = argTr[Up(0)] [l^ . 
For the geometric phase, there are different definitions 
in the literature [lit flTl ITsj . The choice made here 
presents a clear advantage for our purposes, as will be- 
come clear in the following. In order to calculate the 
geometric phase, one should purify the mixed state ap- 
pearing in the Bloch ball. This is done, geometrically, by 
prolongation of the density matrix vector until it reaches 
the unitary Bloch sphere. As a consequence, one gets 
two pure states, pointing in opposite directions (see Fig. 
©)• The geometric phase is given by the weighted sum 
of each pure state's geometric phase. The weights ap- 
pearing in this sum can be obtained from the density 
matrix as follows: one should diagonalize JSJ and express 
it in terms of the two weighted orthogonal components, 
i. e., p = c m \m) (m\ + c n \n) (n|, with (m|n) = 0. Coef- 
ficients c m and c m are the weights of the contribution of 
each state \m) , |n) to the geometric phase (see Fig. 
One possible interpretation of this writing of the density 
matrix is to suppose it comes from the trace of an entan- 
gled state written in the Schmidt's decomposition with 
respect to one of the qubits . 

As a general rule, the total phase for mixed states 
is still a combination of a dynamical and a geometrical 
phase. As will be seen in more details in the following of 
this paper (Section IV), taking as an example a 2tt cyclic 
evolution around a fixed axis, the sum of the geometrical 
and dynamical phases always leads to a total phase gain 
of 7r, as in the pure state case. 

Rotations around one fixed axis studied above can lead 
to some ambiguous conclusions, as pointed out in [l2^ : 
is the effect of ir dephasing a property of the two di- 
mensional rotation group SU(2) or a property of gen- 
eral three dimensional rotations performed by the group 
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FIG. 2: Bloch ball representation of an one qubit density 
matrix and purification procedure. 




FIG. 3: geometric picture of the weighted sum of the tra- 
jectories corresponding to the extension of the density matrix 
vector to two pure states lying on the unitary Bloch ball. The 
areas considered are the marked ones. Coefficients c m and c n 
are the weights of each dashed area. 



SO(3)? Would the results of the experiment change if the 
evolution operator acting on the qubits evolves in time? 
This corresponds to changing the magnetic field direction 
during the evolution in such a way that one always have 
a cyclic evolution at the end. In order to answer these 
questions is to make use of MES. This was done in [T^ 
where a particular geometric picture of the MES Hilbcrt 
space is used. More recently, some extensions to this 
study were made, always with the help of MES . In 
[l2^ we made use of the geometric tools developed in p| 
to explicit the origin of the phase obtained by a MES in a 
cyclic evolution. An interesting property of MES is that, 
as will become clear in the following, it does not gain, in 
a cyclic evolution, neither a geometric nor a dynamical 
phase. However, one can still observe a it dephasing for 
some evolutions, identified as being from topological ori- 
gin. In the following, a more complete geometric picture 
of the space of two entangled qubits is presented, and it 
is shown how different types of phases can appear for non 
MES. 
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III. GEOMETRIC REPRESENTATION OF TWO 
QUBITS 

In this section, we summarize the results presented in 
[l3j where a geometric representation of the space of pure 
two qubit states in presented. This is done with the help 
of Hopf fib-rations, as in the case of one qubit states. The 
main idea is to find a way which is analogous (or at least 
connected) to the Bloch sphere representation presented 
above. 

A. Two Qubits Representation 

Consider a general pure entangled state with arbitrary 
coefficients a, /?, 7 and 5: 

|Vo) =a|00}+/3|01)+7|10} + <y|ll>. (7) 

State (0 is normalized, and can be represented by a S 7 
sphere in R s . A tempting but insufficient way to repre- 
sent (|7J| geometrically is to make a partial trace of one 
of the qubits and study the properties of the resulting 
density matrix. This density matrix can be, as seen be- 
fore, represented in the Bloch ball. This naive method, 
even if dependent on the degree of entanglement of the 
initial state, lacks information: the resulting mixed den- 
sity matrix could originate from an infinite number of 
entangled states with the same degree of entanglement. 
However, with the help of a map defining a Hopf fibra- 
tion, one can combine to this reduced density matrix an- 
other space where the rest of the needed information can 
be found. This other space is the fiber space, defined 
as S 3 /Z2 = 50(3). The fiber space can be represented 
as another S 3 sphere of radius ir with opposite points 
(points differing by a global phase) identified. The com- 
bination of the Bloch ball and the fiber completes the 
Hopf map, where S 7 is described by a basis 5 4 and a 

fiber S 3 via the relation S 7 ^ 5 4 (see Fig Q). It is the 
5 4 sphere's coordinates (base) which are related to the 
Bloch ball. They are given by 

X = (cr a )* = 2Re(a* 7 + (3*5) 

Y = (<r„)* = 2Im(a*7 + 0*6) 

Z = (a z )* = \af + \[3\^\5{ 2 ~\ 1 f (8) 

C r = 2Re(a<5 - £7) 

d = 2lm(aS-p-f). 

The first three coordinates can be identified to the Bloch 
ball coordinates. The last two also appear in a way 
in the Bloch ball: they are connected to its radius as 
r = y/l — C% — Cf . Moreover, C r and C< are related to 
the concurrence C defined by Wooters [2(j by the relation 
C = \C r + id\ = 2\aS — (3j\. The concurrence is a mea- 
surement of the degree of entanglement for this system. 
It assumes the value of 1 for MES and for separable 
states. 



Some particularities of this representation that can be 
directly verified are worth being mentioned. For MES 
(C = I), the Bloch ball reduces to a point. MES are 
thus completely characterized by SO (3). This result is 
related to the fact that all maximally entangled states are 
connected by a local rotation. The reduction of MES to 
50(3) is also the reason why one cannot have geometric 
phases for MES. The global phase in this case should 
come only from the properties of SO (3). For separable 
states (C = 0), the Bloch ball is the Bloch sphere, and 
the same happens to the fiber, so that one ends up with 
two independent Bloch spheres, one for each qubit, as 
expected. In the intermediate case, i.e., for states with 
< C < 1, as mentioned before, one needs 50(3) and the 
Bloch ball to completely represent the entangled state. In 
this case, to each point in the Bloch ball is associated a 
fiber 50(3). 

Now we should gain some insight on how arbitrary en- 
tangled states are represented in SO (3). This can be 
done with the help of Fig. (J3J, where 5*0(3) is repre- 
sented as a sphere of radius 7r with opposite points iden- 
tified. Each point of the 5*0(3) sphere appearing in Fig 
|[SJ is an entangled state. Choosing the origin to be state 
|M/(0)) = a 1 00) +6 1 11) , it is our initial state. Applying a 
rotation to one of the qubits about a given axis displaces 
it in the sphere in the following way: the state performs 
a trajectory in the direction of the applied rotation. The 
angle of rotation is the distance of the final state to the 
the center of the sphere (initial state). In this way, a ro- 
tation about the z axis performed in the first qubit, for 
example, creates a displacement along the z axis of the 
sphere, leading to state R(z, 9) ® I 1*(0)). When 6 = n, 
the state produced can be represented in the north pole 
of the sphere, i.e., ia |00) —iS\ll). Analogously, by ap- 
plying such local rotations, all the other entangled states 
with the same degree of entanglement can be generated 
and represented, as depicted in Fig. JSJ. 

For the purposes of this paper, the geometric represen- 
tation described above is quite useful. We are interested 
in state evolutions that do not change the degree of en- 
tanglement of the initial state and are applied to one 
of the qubits. Such evolutions can thus be geometrically 
represented in two equivalent ways. Defining, in a general 
pure two qubit system ^ ■ \i1j2}, the first and second 
qubits by the subscripts 1,2, let's consider that the evo- 
lution operator acts on the first qubit only. By tracing 
with respect to the second (fixed) qubit, the density ma- 
trix vector relative to the first (evolved) qubit describes 
a trajectory in the Bloch ball. This trajectory leads to a 
geometric and a dynamical phase. At the same time, one 
cannot gain information from the fiber space since, as 
noted above, to each point in the Bloch ball corresponds 
a different fiber space. Alternatively, one can choose to 
make a trace with respect to the evolved qubit. In this 
case, the density matrix vector of the reduced density 
matrix of the fixed qubit doesn't move in the Bloch ball. 
There is no geometric nor dynamic phases. However, 
the fiber space is always the same, and the entangled 
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state performs an evolution in this space. This evolution 
leads to a total global phase of topological nature that de- 
pends on the trajectory described by the entangled state 
on SO(3). In the next section, we study the details of 
these different phase dynamics during an entangled state 
evolution. In order to do so, we first set notations to 
see more clearly the reduced density matrix vector state 
evolut' ' " 

SO(3) 




C=0 0<C<1 C=1 

S 2 xS' S 2 xSO(3) 50(3) 



FIG. 4: Geometric representation of a pure two qubit sys- 
tem as a Bloch ball obtained from the trace with respect 
to one of the qubits and a sphere representing SO(3). This 
representation is a function of the concurrence C of the sys- 
tem. In the limit cases of MES (C = 1) and product states 
(C = 0), this representation reduces to 5*0(3) or a product of 
two Bloch spheres, respectively. For a general non MES state 
(0 > C > 1), the pure qubit state is represented by the Bloch 
ball and SO (3) 




FIG. 5: A possible representation of SO(3): a ball of radius 
7r with opposite points identified. Entangled states are points 
in this ball and displacements on it correspond to rotations 
applied to the entangled states. States in the border of the 
sphere are obtained by rotation of one of the qubits of the 
initial state (placed in the center of the sphere) of ir. 



B. The Reduced Density Matrix in the Bloch Ball 

Let's take as an example state 

H'e) = v^cos - |0i0 2 ) - VATsin - |0il 2 ) (9) 



+ V^sin- |li0 2 ) + V^icos - |lil 2 ) . 

The subscripts define the first and the second qubit, and 
they will be omitted in the following. The concurrence 
of © is 2VA G Ai. The angle 9 has here a geometric inter- 
pretation: coordinates X, Y and Z of the state vector in 
the Bloch ball are given, using the results of the previous 
section, by 

X = (1 - 2A )sin(9 

Y = (10) 

Z = (2A Q - 1) cos 6>. 

The radius of the Bloch ball is r = |2A D — 1| and all states 
lying on its surface have the same degree of purity (and 
come from two qubit states having the same amount of 
entanglement). The angle 9 is the angle the density ma- 
trix vector, represented in Fig. lf2"|. makes with the z 
axis. In order to gain some geometric insight, we can 
start by writing the general reduced density matrix and 
studying some specific cases. The reduced density ma- 
trix, after tracing with respect to the second qubit, is 
given by 

pe = (A o cos 2 ^ + A lS in 2 ^)|0) (0| + 

(Ai cos 3 | + Ao sin 2 |)|1><1| (11) 
(Ao-Ai)sin0(|l} <0| + |0)(1|. 

For 9 = 0, we have the entangled state |^o) = 
•v/A^|00) + vi — A |11). The one qubit density ma- 
trix corresponding to this state is po = A Q |0) (0| + 
(1 — A ) |1) (1|. The density matrix vector is thus par- 
allel to the z axis. If 6 = ir/2, we have |^ , 7 r/ 2 ) = 

^(v%|oo) + VT - X|oi) - VK\w) + a/T - — \n)). 

The reduced density matrix in this case is p^li = 
1/2 (|0) (0| + |1) (1|) + ^ii (|0) (1| + |1) (01), and the 
vector representing it in the Bloch ball lies in the equa- 
tor. We can see that writing JJjJ as a function of 9 makes 
the geometric representation of the Bloch vector more di- 
rect. Let's consider now rotations performed in the first 
qubit of the entangled state. Such rotations also rotate 
the reduced density matrix vector exactly as it happens 
to pure states. Taking once again as the initial state 
I*) = VK\00) + y/l - A |11) {9 = 0), its density matrix 
vector lies in the z direction. Rotations of the first qubit 
of an angle 9/2 with respect to the x axis via the transfor- 
mation R{9/2) lead to a vector making an angle 
9 with the z axis. Rotations around other directions in 
space can be done in an analogous way. 

IV. CYCLIC PHASES FOR PURE TWO QUBIT 
STATES 

We can now investigate the possible origins of the 
phase gained by a general pure two-qubit entangled state 
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after a cyclic evolution. We again choose, for simplicity, 
to consider evolution operators acting on one qubit only, 
for example, the first one, as defined in ©. After a time 
interval T the final state is brought to the initial one, 
apart from a phase factor. The initial state of the system 
is of the form of (0. This state, as seen in the previous 
section, can be represented geometrically with the help 
of a Bloch ball of radius proportional to C and of SO(i). 
The general matrix form of the evolution operator is: 

jy _ e ~wt _ ( cos | — sm | —in- sin | \ 
\ —in+ sin | cos | + in z sin | J ' 

(12) 

where H = ha, with a = a x x + a y y + a z z. <7, are the 
Pauli matrices, h is the unitary vector with coordinates 
h = n x x + n y y + n z z and h± = h x ± ih y . Eq. I|12(l 
above represents a rotation around an axis h in space. 
The evolution operator described in (|12fl acts only on 
one qubit, and the total evolution operator acting on 
both qubits is given by the product U ® I or I (g) U. The 
total evolution is given by a sequence of operators of the 
type (|12fl acting on one qubit of the initial state Q , each 
one of them for a fixed time interval t. 

In the following subsections, we study the particular 
case of MES and then generalize our results to states with 
a variable degree of entanglement. The case of MES has 
already been studied in [l2] , and is of particular interest 
because only in this case the origin of the cyclic phase is 
purely topological. For non MES, we will see that geo- 
metric and dynamical effects can harmonize themselves 
so as to create a cyclic phase of wr, equal to the topo- 
logical phase. 



A. Maximally entangled states 

MES, as stated before, are a very particular case in 
what concerns the global cyclic phase. Since the degree 
of entanglement is maximal, its reduced Bloch ball is of 
null radius. It is completely represented by SO(3). For 
this reason, there is no geometric phase. We can also 
easily check that, using the definition for the dynamical 
phase 4>d 0> there is no dynamical phase associated to 
a MES in a cyclic evolution if one uses a sequence of op- 
erators of the form l)12H . Nevertheless, depending on the 
choice of the set of operators applied to the initial state, 
one can either have a cyclic evolution leading to no global 
phase change ("plus" trajectory) or another one where 
a global phase of n appears ("minus" trajectory). The 
existence of these two types of trajectories of MES in 
SO(3) comes from the biconnectcd nature of this group 
[Tflj . The "plus" and "minus" trajectories belong to two 
different homotopy classes. The difference between both 
classes of trajectories can be seen geometrically from the 
spherical representation of SO (3): if one places the ini- 
tial state at the origin of the SO(3) sphere, the border 
of this sphere represents the space of states orthogonal 
to the initial one. They are mathematically represented 



by U\^(0)} = R(ir,n) <g> I|*(0)). Whenever the initial 
state reaches this space, there is a discontinuity in the to- 
tal phase. If the evolved state crosses this space an odd 
number of times, a total global phase of tt appears after 
the cyclic evolution, defining the "minus" type trajecto- 
ries. In the "plus" type, it crosses the same space an even 
number of times. The parity of the number of times the 
evolved states crosses the space of states orthogonal to 
the initial one is responsible for the topological difference 
between trajectories and for the different phase gain. An 
example of a "plus" and a "minus" trajectory was stud- 
ied in 12] . They correspond to two distinct sequences 
of rotations applied to one of the qubits. The "plus" 
trajectory is performed by a sequence of operators of the 
form (|12|l generating the trajectory represented by the se- 
quence A-^B^F^D^A and the "minus" by the 
sequence of operators generating A^B^F^E^A, 
where A = —A and E = —E. A, B, D, E, F are points 
in SO(3) and both trajectories are depicted in Fig.©. 
Recall that rotations in SO (3) are represented as dis- 
placements in the direction of h of a distance given by 
the angle of rotation 6. Notice that, for the minus type of 
trajectory, represented pictorially by the solid line, after 
crossing the border of the sphere in F, the state reappears 
in the opposite side of the sphere and continues its trajec- 
tory from F. For each part of the trajectories mentioned 
above, the rotation operator acts for a time t = 2tt/3 and 
the only difference between the evolution operators 1|12|) 
performing each part of the trajectory is the direction to 
which h points. These orientations (n x ,n y ,n z ) read: 



A—>B 


v/T73(-i,-i,-i) 


B -> F 


073(1,-1,-1) 


F^D 


v/l/3(-l,-l,l) 


D^A 


Vl/3(-l,l,l) 


F -> E 


^/l/3(-l,-l,-l) 


E^A 


v/l/^l,-!,-!) 



As can be seen, the axis of rotation varies in each one of 
the trajectories, lifting the ambiguity on the origin of the 
7r phase: it is manifestly a property of SO (3). The case 
of MES subjected to a sequence of unitary evolutions 
of the type of Eq. I|12|) is an example where there is 
no geometric nor dynamical phase. The only possible 
origin of the global phase gain is thus topological. We 
see now how these results are affected by the degree of 
entanglement. 

B. States with an arbitrary degree of entanglement 

We study now the most general case of two-qubits with 
an arbitrary degree of entanglement. The limit case of 
a product state (that should reproduce the results ob- 
tained for two independent qubits precessing in their 
Bloch sphere) appears naturally from the general results. 
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FIG. 6: Two possible trajectories belonging to two different 
homotopy groups in 50(3): sequence A-^B^F^D^A 
is a plus type trajectory that does not cross the space of states 
orthogonal to the initial one (A). Sequence A — > B — » F — > 
E — > yl is a minus type trajectory, leading to a global phase 
gain of 7r. Note that, in this geometric view, points A and 
A are both in the center of the sphere, since this point is 
opposite to itself. 

Consider an arbitrary initial state of the form of Q 
where the first qubit only evolves, under the action of 
a given Hamiltonian. In principle, one should expect, 
when tracing with respect to the second (fixed) qubit, 
to have a geometric phase, since the resulting Bloch ball 
corresponding to the reduced density matrix has a non 
null radius and the qubit describes a trajectory on it. 
This trajectory leads to the appearance of the expected 
geometric phase. At the same time, there should be a 
dynamical phase, and some natural questions are: what 
is the total phase of the non MES after one cyclic evo- 
lution? What are the contributions of the dynamical, 
geometrical and topological phases to this total phase? 
How does it depend on entanglement? A way to start 
to answer this question is to apply to the non MES the 
two sequences of transformations that were performed to 
the MES in the previous section and calculate the total 
phase with the help of Taking as initial state, for 
example, 

|*(0)> = v / ^|00) + v / l-A o |ll) (13) 

there are two equivalent scenarios, that must lead to the 
same final result since they depend on the pure entan- 
gled state only: tracing (|13[) with respect to the evolved 
or the fixed qubit. If one traces out the evolved qubit, all 
the phases should come from the trajectory of the non 
MES in SO(3). By tracing with respect to the second 
(fixed) qubit, the phase information must come from the 
reduced density matrix vector's trajectory in the Bloch 
ball and the dynamical phase. Let's consider first the 
case where one traces out the evolved qubit. Defining 
(|T3|> as the center of the sphere representing SO (3), we 
have that, after applying the same sequence of rotations 
discussed in the previous subsection, there are still two 
types of trajectories (plus and minus) in 50(3) that are 
represented exactly in the same way as Fig.©. In gen- 



eral, we therefore have that the total phase after the cyclic 
evolution is the same as for the MES case regardless of 
the value of A Q . 
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FIG. 7: Total phase for trajectories "minus" (above) and 
"plus" (below) for three values of the parameter A : 0.3 (dot- 
dashed line), 0.4 (dashed line) and 0.48 (continuous line). The 
curves become more sharp and tend to discontinuity as the 
value of A increases. For A = 0, 1 (product states), it is 

cri^ron Yi\t a cfmitrnt lino ac pvnpftufl 




FIG. 8: Geometric phase for a cyclic evolution around a fixed 
axis as a function of the angle 6 and of parameter A , propor- 
tional to the degree of entanglement. 

Indeed, the existence of the two types of trajectory 
doesn't depend on the degree of entanglement. The to- 
tal phase can be calculated from (J). Fig. J7J) shows 
it for three possible values of A„ = 0.3, 0.4 and 0.48 for 
both types of trajectory. We see that for the three values 
of A one can, effectively, define two types of trajectory. 
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However, the closer we are to MES (A D = 0.5), the more 
abrupt the phase gain becomes, until it reaches a discon- 
tinuity for A G = 0.5. In this limit case, the phase is either 
zero or 7r during the whole trajectory and the discontinu- 
ity represents the gain of the topological phase. As seen 
before, it is a consequence of state's \^(t)) crossing the 
border of the 5*0(3) sphere in the middle of the trajec- 
tory (t = 47r/3). FigQ shows that the same happens for 
non MES. 

Let's consider now the same trajectories but trace with 
respect to the fixed qubit. This results on the reduced 
density matrix of the evolved qubit. This density matrix 
performs a trajectory in the Bloch ball. As mentioned 
above, the existence of the two types of trajectory is a 
mathematical result, valid for the entangled pure state, 
that must be reproduced in the two situations: tracing 
with respect to the evolved or the fixed qubit. By tracing 
with respect to the fixed qubit, one can calculate the 
global phase from its definition (Q. In order to do so, we 
should first calculate the scalar product: 

SP = cos ^ + i(Xn x (t) + Yn y (t) + Zn z (t)) sin ^. (14) 

In the expression above, X, Y and Z are the coordinates 
of the Bloch vector of the reduced density matrix and 
the vector n was taken as a function of time, as proposed 
in [Tflj . In the example studied, n(t) suffers discontinu- 
ous changes corresponding to the different axes of rota- 
tion. We see that Eq. l|14|) is completely determined by 
the coordinates of the reduced density matrix vector in 
the Bloch ball and the same result could be obtained by 
calculating the mixed state total phase = argTr[[//5]. 
One can see that for the particular initial state (|13fl and 
the sequence of rotations under consideration, the scalar 
product never reaches the value of zero, so there is no 
discontinuity in the phase. This result has already been 
seen in Fig.Q as being a consequence of the trajectory 
of the entangled state in 50(3). However, Eq. l|14|l shows 
clearly that the same result can be obtained from in- 
formation on the Hamiltonian and the reduced density 
matrix only. Ecr. (|14|) also shows that a total ?i7r phase 
still appears. This result can be seen as being a conse- 
quence of the parity of the number of times the initial 
state crosses the border of the 50(3) sphere or of the 
area of the trajectory of the reduced density matrix in 
the Bloch ball combined to the dynamical phase. 

We can now calculate the contribution coming from 
the dynamical and the geometric phases of the reduced 
density matrix. The calculation of the dynamical phase 
is performed from its definition It is given by the 
integral 

1 f* 

- J (n x (t)X + n y (t)Y + n z (t)Z)dt, (15) 

that can be seen as the integral of the scalar product 
between the density matrix vector and the direction of 
rotation hit). For this specific example, the dynamical 
phase is still zero. It is thus the geometric phase that 



explains the existence of the two types of trajectories. 
The two different sequences of evolution operators lead 
to different trajectories in the Bloch ball that are inde- 
pendent of the degree of entanglement. They are shown 
in Fig.©. The area of each one of this trajectories can be 
calculated, and by weighting them using the techniques 
of Section II, one finds the two possible global phases, 
or 




FIG. 9: Trajectories in the Bloch ball for the two possible sets 
of rotations (a) A — > B — > F — > E — » A (minus trajectory), 
and (b)A^B^F^D^A (plus trajectory). 

We see that either one can obtain the total global phase 
of a non MES entangled state from the topological prop- 
erties of 5*0(3), or by geometric properties of a reduced 
density matrix. This means that by looking at a subsys- 
tem of the entangled state gives us information about its 
global phase. 

As a general rule, the scalar product (|14f> and the dy- 
namical phase (|15|) can tell us all about the types of 
phases appearing in an evolution. To illustrate that, 
let's study a simpler Hamiltonian, describing a rota- 
tion around only one fixed axis, say, the z axis. For 
a general state of the type ©, l|14|) depends on 9 as 
SP = cos | + i(2A — 1) cos 9 sin 4. The dynamic phase 
for this case is 4>d — 7r(2A — 1) cos 6*. Notice that, apart 
from a scaling factor depending on the degree of entan- 
glement (2A — 1), it has the same form as for pure states. 
The cyclic geometric phase also depends on 9 and A . It 
is given by 4> g = ir + n(l — 2A„) cos 9, plotted in Fig. (JHJ. 
The connection of this phase to the pure state case is not 
direct. However, some limit cases are reproduced. First, 
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notice that for A D = 1/2 (MES), 4> g and <fid vanish, as 
expected. Also, for separable states (A Q = 0,1), <p g os- 
cillates with 9 as in the pure state example discussed in 
Section II. At the same time, in this limit, the total phase 
is a linear function of time, as one should expect. The 
dependence of the phases on 9 of 4>g is also interesting: in 
the limit 9 = tt/2, <fid and 4> g show that the total global 
phase is of geometric nature only, a result analogous to 
the one found for one qubit states. However, for 9 = 0., 
the analogy with the one qubit case breaks down. The 
geometric phase is given by (f> g = 2tt\i in contrast to 
zero for the case of one qubit. Nevertheless, there is no 
contradiction, since in the discussed case of mixed states 
the dynamical phase doesn't account for the total phase. 
It assumes the value of ir(2\ — 1), so both add together 
to 7r, as expected. It is also easy to verify that for an ar- 
bitrary value of 9, <p g + (f>d = 7T, while the proportions of 
geometrical and dynamical phases are determined by the 
angle 9. By looking now again at the case where the first 
qubit is traced away, for an evolution representing a rota- 
tion around one fixed axis in space, all trajectories cross 
the border of the SO(3) sphere, accounting for a global 
phase of 7T. This is why we must observe <f> g + 4>d = 7r '. 

The discussion above generalizes the results of [TjI and 
gives a geometric interpretation of the nature of the cyclic 
phases for entangled two-qubit states with an arbitrary 
degree of entanglement. This is done with the help of 
geometric phases represented in the Bloch ball for den- 
sity matrices. As for the case of MES, the results pre- 
sented can be, in principle, tested experimentally To 
do so, one needs to throw the non MES in an interfer- 
ometer, so that the global phase can become a relative 
phase. A possible wa y to do that is using photon interfer- 
ence, as proposed in [13. Recently, other proposals have 
appeared aiming to show the same effects in a system 
of entangled spin pairs 0. Both proposals need slight 
modifications in order to account for the phase effects 



described in this paper. This modification concerns only 
the initial state, that now can have an arbitrary degree of 
entanglement. Taking the example of polarized twin pho- 
tons, it has already been demonstrated that polarisation 
entangled photon pairs with an arbitrary degree of en- 
tanglement can be generated [2l|. As shown above, the 
same set of transformations performed in MES lead to 
the two different types of trajectories also for non MES. 
This means that the same set of wave plates used in 
can be employed to test the predicted results. 



V. CONCLUSIONS 

To conclude, the different types of phases appearing 
in a cyclic evolution of a pure bipartite entangled state 
were studied. This is done from a geometric point of 
view. It is shown that the interplay between topologi- 
cal, dynamical and geometrical phases always lead to a 
global phase of nir. The parity of n, can have different 
origins, depending on how one choses to geometrically 
represent the pure entangled state: by considering evolu- 
tion operators acting on one qubit only and tracing with 
respect to the fixed qubit, trajectories of the type plus 
and minus are a consequence of geometric and dynamical 
phases of the evolved qubit. Alternatively, still supposing 
that only one qubit only evolves and tracing it out, the 
parity of the total phase is a consequence of the number 
of times an initial state, placed at the origin of SO(3), 
crosses the border of this sphere. As a consequence, it 
is shown that mixed states coming from the trace of a 
pure entangled state with respect to one of the particles 
may play a crucial role in describing the total phase gain 
of the pure entangled states itself. The author wishes to 
thank R. Mosseri for inspiring discussions and I. Lorgere 
for a critical reading of the manuscript. 
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